Dynamical conductance in the two-channel Kondo regime of a double dot system 
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We study finite-frequency transport properties of the double-dot system recently constructed to 
observe the two-channel Kondo effect [R. M. Potok et al, Nature 446, 167 (2007)]. We derive an 
analytical expression for the frequency-dependent linear conductance of this device in the Kondo 
regime. We show how the features characteristic of the 2-channel Kondo quantum critical point 
emerge in this quantity, which we compute using the results of conformal field theory as well as nu- 
merical renormalization group methods. We determine the universal cross-over functions describing 
non-Fermi liquid vs. Fermi liquid cross-overs and also investigate the effects of a finite magnetic 
field. 
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I. INTRODUCTION 

In the past few years semiconducting quantum dots 
have been in the focus of intense research. This research 
has mostly been motivated by their possible application 
in future microelectronics: these devices behave as tun- 
able artificial atoms attached to electrodes,— they can 
be used as single electron transistors,^ furthermore they 
also serve as a playground to model and study (artifi- 
cial) molecular transport in a very controlled way. They 
display a number of correlation-induced effects like the 
Coulomb blockade or the Kondo effect^-'* and they can 
also give rise to exotic strongly correlated states.^ Very 
importantly, quantum dots can also be used to build 
quantum bits with the electron spin providing the neces- 
sary degree of freedom for quantum-computation.™ 

Nevertheless, maybe the most fascinating application 
of quantum dots is their possible use to realize quan- 
tum phase transitions between different correlated states. 
Several such transitions have been proposed: under spe- 
cial circumstances the transition between the triplet and 
the singlet state of a dot can be a true quantum phase 
transitioUfii^ although in most cases this transition be- 
comes just a cross-over.^ Dissipation can also lead to 
a quantum phase transition where the charge degrees 
of freedom of the dot become localizedi ^-'^d^ Unfortu- 
nately, these phase transitions have a Kosterlitz-Thouless 
structure and are - in a sense - 'trivial' quantum phase 
transitions. Using multi-dot systems, however, it is 
also possible to realize generic quantum phase transi- 
tions, where the transition point represents a true crit- 
ical state characterized by anomalous dimensions and 
a singular behavior. These critical states are generic 
non-Fermi liquid states in the sense that they cannot 
be described in terms of conduction electron quasiparti- 
cles even at the Fermi energy.^* The prototypes of these 
generic quantum impurity states are the two-channel 
Kondo modeii^ and the two-impurity Kondo model. 



Some years ago Matveev proposed that the two-channel 
Kondo model could be realized by charge fluctuations at 
the charge-degeneracy point of a quantum dot.— How- 
ever, Matveev's mapping assumes a vanishing level spac- 
ing and with present-day technology it has been impos- 
sible to reach this state so far. However, a few years ago 
Oreg and Goldhaber-Gordon proposed to realize the two- 
channel Kondo state through a double dot systemfi^ and 
after several years of work this two-channel Kondo state 
has indeed been observed in a pioneering double dot ex- 
periment at Stanford. For the realization of the other 
prototypical non-Fermi liquid state, the two-impurity 
Kondo state, a somewhat similar multi-dot setup has 
been proposed recently^22 

Fig. [1] shows the double dot device suggested by Oreg 




FIG. 1: Two-dot device: the small dot in the center couples 
to a large dot (2) and to a left and a right lead {IL and IR) 
via the hopping amplitudes: vl and vr. The small dot has 
a large level spacing, and the large dot is characterized by a 
vanishing level spacing, while both dots are in the Coulomb 
blockade regime. As a result, only spin exchange is possible 
between the dots. 
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and Goldhaber-Gordon, which has since been used to ex- 
perimentally reach the two-channel Kondo fixed pointi^. 
This set-up consist of a small dot coupled to a large dot 
(2) and two leads (li and IR). The small dot is tuned 
to the regime where charge fluctuations are suppressed 
and it has only one extra electron on it. The level spac- 
ing, des, of the small dot and its charging energy ~ Ecs 
are assumed to be much larger than the temperature, 
Scs , Ecs ^ T, so that below the scale D charge fluctua- 
tions on the small dot are suppressed and the only role 
of this dot is to provide a spin. The size of the large 
dot, on the other hand, is chosen in such a way that its 
charging energy and level spacing satisfy Ec2 > T > 6e2. 
This implies that this dot is also in the Coulomb blockade 
regime while the electronic states on it form a continuum 
of electron-hole excitations. Therefore, electrons on the 
large dot form a bath that can exchange spin with the 
small dot while electrons cannot jump out of it^^ as it 
is also indicated in Fig. [T] In the limit of small tunnel- 
ing amplitudes, apart from some irrelevant and potential 
scattering terms, this double dot system is described by 
the following simple two-channel Kondo Hamiltonian, 



(1) 



The operator ip2 describes electrons on the large dot. In 
the contimmm limit, Se2 — + 0, it is defined as 
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y ao-(e) de 



(2) 



with a^{e) the annihilation operator of a conduction elec- 
tron of energy e and spin a on the large dot, satis- 
fying the anticommutation relation: {ao-(e), a^/(e')} = 
Saa' S{e — e'). The operator in Eq. dT]) is a suitably 
chosen linear combination of electrons on the left and 
right lead-electrodes. 



VL^L + VRlpf 



(3) 



with vl and vji the hopping amplitudes between the dot 
and the left and right electrodes, respectively. The left 
and right fields ipL /r are defined similarly to Eq. ([2]) , 



(4) 



with Ci^/ji ,j{e) the annihilation operator of a conduction 
electron of energy e and spin a on the left /right lead. 

We remark that, strictly speaking, the Kondo Hamil- 
tonian above is only accurate in the limit of small tun- 
neling, while in the experiments the tunneling rates were 
quite large in order to boost up the Kondo temperature.-iSi 
Therefore, to study the region far above Tk, an Ander- 
son model-type approach would be needed that also ac- 
counts for charge fiuctuations of the small dotfSZ, Never- 
theless, our Kondo model-based approach captures accu- 
rately the universal cross-over functions in the region of 
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FIG. 2: Top: Sketch of the conductance through the small 
dot divided by its maximum value. Go, as a function of tem- 
perature. For Ji = J2 a ~ VT singularity emerges, while 
for Ji 7^ J2 a Fermi liquid is formed at a scale T* , and the 
conductance crosses over to a very small or a large value, 
with a characteristic Fermi liquid scaling, ~ (T/T*)^. Bot- 
tom: Sketch of the "phase diagram" of the two-channel Kondo 
model. 



interest, i.e. around and far below the Kondo temper- 
ature, provided that both dots are close to the middle 
of the Coulomb blockade regime. To account for devi- 
ations from the middle of the Coulomb blockade valley, 
one could break the particle-hole symmetry of Eq. p]) 
and add potential scattering terms to it. 

The quantum critical state arises as a result of the 
competition of channels 1 and 2 to form a singlet with 
the dot spin S. Depending on the values of the dimen- 
sionless couplings, Ji^2, two situations can occu n^^i^^ (a) 
For Ji < J2 the spin of the small dot forms a Kondo 
singlet with electrons on the large dot that screen the 
spin at an energy scale T* . In this case, to promote a 
conduction electron between the left and right leads one 
needs to break up the Kondo singlet and pay an energy 
T* , and therefore transport through the small dot is sup- 
pressed at low temperaturesi ^^'^^ (b) For Ji > J2, on the 
other hand, the spin of the small dot is screened by elec- 
trons in the leads. This correlated singlet state serves as 
a 'bridge' and helps the lead electrons to propagate be- 
tween the left and the right sides with a small refiection 
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probability and is thus characterized by a conductance of 
the order of the quantum conductance, 2e^ /h. In both 
cases a Fermi hquid state is formed below the scale T*, 
which can be characterized by simple phase shifts at the 
Fermi energyi^ 

Interestingly, for Ji , J2 J the scale T* vanishes as 
T* ~ (Ji — J2)^, and a non- Fermi liquid state emerges 
below the Kondo scale, Tk ~ D e^^l \ with the cut-off 
defined as D = min{(5es, Ecs^Eci)^^ This so-called two- 
channel Kondo state is characterized by a conductance 
that is about half of the quantum conductance at very 
low temperatures, and has a ~ yjT /Tk singularity for 
T <C TV-^^ This state is - in a sense - a quantum critical 
state: although it is just a single point in the parameter 
space, it separates two stable Fermi-liquid phases, and 
it influences the behavior of the double dot system over 
the whole regime, T* < T,uj < Tk for Ji ^ J2. How- 
ever, as we shall see later, the scaling properties usually 
associated with the two-channel Kondo fixed point it- 
self are restricted to a somewhat smaller energy range, 
^JT*Tk < T,uj < Tk ■ The characteristic features of the 
temperature-dependence of the DC conductance and the 
schematic phase diagram are sketched in Fig. [2l 

The purpose of the present paper is to investigate dy- 
namical transport properties of the above set-up and de- 
termine how the two-channel Kondo behavior and the 
presence of a quantum critical point at Ji = J2 manifests 
itself in the AC conductance through the dot. For this 
purpose, we shall derive an expression for the AC conduc- 
tance in the linear response regime that relates the con- 
ductance to the so-called composite Fermions' propaga- 
tor at any temperature and frequency. Making use of this 
simple formula, we shall evaluate the AC conductance 
for frequencies T <^ u using numerical renormalization 
group methods. We shall also determine the universal 
cross-over functions that describe the non-Fermi liquid 
vs. Fermi liquid cross-over for T* <C Tk- As we show, 
the AC conductance exhibits features that are qualita- 
tively similar to the finite temperature DC conductance, 
sketched in Fig. [51 In addition, we shall also investigate, 
what conclusions we can draw regarding AC properties 
based upon the predictions of conformal field theory, and 
use this approach to obtain the universal scaling of the 
conductance in the regime T* <ti u^T <^ Tk ■ 

The paper is organized as follows. Section [III provides 
the details of the derivation of the AC conductance for- 
mula for the two-channel Kondo model. In Section Hill we 
present some analytical considerations based on confor- 
mal field theory concerning the universal scaling proper- 
ties of the linear conductance and of the eigenvalue of the 
so-called on-shell T-matrix. Section HVl comprises our nu- 
merical renormalization group results for the composite 
Fermions' spectral function and the linear conductance 
in case of channel anisotropy and in the presence of a 
magnetic field. At last our conclusions are summarized. 



II. KUBO FORMULA AND COMPOSITE 
FERMIONS 

Let us start our analysis with the derivation of a sim- 
ple expression for the AC conductance in terms of the 
so-called composite Fermion operators/^ For this pur- 
pose, we first couple an external voltage to the dot and 
introduce a time-dependent chemical potential difference 
between the left and right electrodes: 



Hv = V{t) Q = e V{t) {Nr - Nl) 



(5) 



with Nr and the number of electrons in the left and 
right leads, respectively. 



L/R 



The current operator can be defined as the time deriva- 
tive of Q, I{t) = i [H,Q{t)] = i [Hint:Qit)]. This com- 
mutator is easily shown to give 



+ h.c. 



(6) 



where V'l denotes the decoupled electron field of the leads 



(7) 



and we have introduced the so-called composite Fermion 
operator. 



(8) 



The operator F'^ has spin 1/2 and charge 1, and it cor- 
responds to the 'universal part' of the electron localized 
on the small dot. 

Close to equilibrium, the current through the dot is 
given by the Kubo formula 



(/(<)) = j G[t~ t') V{t') dt' , 
G{t-t') = t{[m,Q{t')])0{t^t'), (9) 

with G{t — t') the conductance. Differentiating with re- 
spect to time and then taking the Fourier transform we 
obtain the relation 



(10) 



where Qf^ denotes the retarded current-current correla- 
tion function and ^ is a real constant 



(11) 



Thus the real and imaginary parts of the conductance 
are given by 



Re{G(cj)} 



1 



ImiGfiito)} 



(12) 



Im{G(c^)} = {Re{g^,{uj)}-Re{gfj{0)}) . (13) 
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In general, it is not so simple to compute the correla- interacting Fermions. This observation allows us to write 
tion function In our case, however, the field tpi is Quit) as 
completely decoupled from the spin and describes non- 



,2 ,,2 



where C/^, Q^, , and Q'^ denote the usual retarded, advanced, bigger and lesser Keldysh Green's functions. The 
Fourier transform of this expression simplifies considerably if one uses the fact that the field ipi is non-interacting and 
therefore the corresponding Green's functions become in the large band-width limit 



(15) 



with f{Lu) the Fermi function. Taking the real and imaginary parts of the Fourier transform of Eq. p4p we finally 
obtain: 



Re{GH} = / duj'lmiUuj')} [f {uj' + u;) - f {uj' - u)] , 

Im {G(c^)} - ^Y^Jdw'Re {Uiuj')} [f{cu' + lo) + f{J ~ uj) - 2/(c.')] 



(16) 
(17) 



where we introduced the dimensionless eigenvalue ta{i^) 
of the so-called on-shell T-matrix,^^ which describes the 
scattering of electrons of energy w. 



conductance through the dot in the vicinity of the two- 
channel Kondo fixed point. 

III. ANALYTICAL CONSIDERATIONS 



(18) 



and Go denotes the maximum conductance through the 
dot, 



Go 



2e^ 



h {v 



'AY 



(19) 



Thus the real part of the conductance is related to the 
imaginary part of gp„^ which is essentially the spectral 
function of the composite Fermion, gpuiuj). The latter 
can be determined numerically using the numerical renor- 
malization group method. Then the real part, Re{GP^}, 
can be obtained by performing a Hilbert transformation 
numerically, and the imaginary part of the conductance 
can then be calculated from Re{G^^} by simple numer- 
ical integration. Note that Eqs. (fT5)) and P7| provide 
the linear conductance through the dot for any asym- 
metry parameter at any temperature and any frequency. 
They are thus natural extensions of the formula given in 
Ref. [l7, and are the analogues of the formulas obtained 
recently for the Anderson model. 

Eqs. (Uni) and ^ belong to the main results of our 
paper. We shall use these formulas to compute the AC 



Eq. p6|) allows us to make numerous statements based 
on rather general properties of the two-channel Kondo 
fixed point. From an exact theorem of Affleck and 
Ludwig,'^*' e.g. we know that at the two-channel Kondo 
fixed point (i.e. for Ji = J2 and aj,T — > 0) the 5-matrix 
of the conduction electrons identically vanishes. From 
the relation, S{ijj) = 1 + i t{uj) between the dimension- 
less eigenvalue of the 5'-matrix and the T-matrix wc thus 
obtain 



lim t(uj,T)^i ( Jl = J2) 



(20) 



From this, it immediately follows that at the two-channel 
Kondo fixed point the conductance takes half of its max- 
imum value. 



lim G(w,T) = Go/2, (Ji - J2) . (21) 



The results of conformal field theory"^" also enable us to 
compute the finite frequency conductance for Ji = J2 
and Lj,r <C Tk- In this limit the T-matrix is given by 
the expression^ 



5 



3(7rr)^/^A / du 





-1/2 



(1 



-3/2 



(22) 



where F{u) = F(3/2, 3/2, 1; m) is the hypergeometric 
function, and A stands for the amphtude of the leading 
irrelevant operator: 

^ " vfc ■ ^^^^ 

The value of the dimensionless constant 7 depends on 
the precise definition of Tk- Throughout this paper, we 
shall define Tk as the energy at which for Ji — J2 the 
composite Fermion's spectral function drops to half of its 
value, Im t{uj = Tk) = Ini t{uj = 0)/2. Then, comparing 
the numerical results of Section IIVI to the asymptotic 
m ^ T behavior of the conductance we obtain the value 
7 = 0.093 ± 0.001. Clearly, since the omega-dependence 
enters t{Lj) only in the combination w/T, it immediately 
follows that l-Imt(cj,r)/(ATi/2) is a universal function 
of Lo/T (see inset of Fig. 3). 

In Fig. [3] we show the results obtained by numerically 
integrating Eq. (|22p for a few temperatures. It is remark- 
able that curves corresponding to different temperatures 
cross each-other. This feature is a direct consequence of 
the unusual shape of the universal curve shown in the 
inset of Fig. H 

Note that to construct the universal scaling curve one 
needs to rescale the axes with respect to the temperature 
only, and the precise value of the Kondo temperature ap- 



pears only through the prefactor A. The fact that the 
only relevant energy scale is the temperature is charac- 
teristic of quantum critical points. The imaginary part of 
the T-matrix exhibits a ~ V^M^ cusp for T <^ uj <^ Tk, 
and crosses over to a quadratic regime for lo <ti T. Simi- 
lar behavior is observed in the real part of t(w), shown in 
Fig. m This quantity also shows a characteristic ~ ^/uj 
behavior at frequencies Tk ^ lo ^ T, that crosses over 
to a linear regime for lj <^T . 



Using Eqs.dUl), Eqs.lUS]) and dTT]), both the real and 
the imaginary parts of the conductance can be computed 
by numerical integration. The results are plotted in 
Figs. [5] and [6] for various temperatures. Even though, 
at first sight, the results for the conductivity look qual- 
itatively similar to those for the T-matrix, there is an 
important difference: integration with the Fermi func- 
tions apparently eliminated the aforementioned crossing 
of the curves. Similar scaling curves have been computed 
using conformal field theory results for the differential 
conductance of two-channel Kondo scatterers in point 
contacts^i 




K 

FIG. 3: Imaginary part of the eigenvalue of the T-matrix ob- 
tained by numerical integration of Ea. l|22p . The scale of the 
uj axis is set by the amplitude of the leading irrelevant opera- 
tor, A. The inset illustrates how the curves corresponding to 
different temperatures collapse into one universal curve. 




.41 I I I I 

-0.2 0.2 

FIG. 4: Real part of the eigenvalue t{uj) of the T-matrix pre- 
dicted by conformal field theory. The inset shows the collapse 
to a single scaling curve (obvious from the integral definition) . 
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In the limit, Tk ^ w ^ T, the conformal field theory 
also predicts that the AC conductance scales as 



Re G{lu) 



Go 



1 — a 



ImG(w) « ^ a sig-n{uj) J ]^ , 




(24) 



with a — 2.53 ± 0.06 a universal constant of order of 
unity. The fact that the coefficients in the two equations 
above are both equal to a, follows from the observation 
that G{uj) is analytical in the upper half-plane. 

For Ji J2 a new Fermi liquid scale, T* (mentioned 
earlier), emerges, but one can still make many statements 
based upon the fact that the leading relevant and irrel- 
evant operators have scaling dimensions y+ = 1/2 and 
?/_ = —1/2, respectivelyjAS. As a consequence, in the 
vicinity of the two-channel Kondo fixed point (T* <C Tr-) 
the conductance becomes a function of the form: 



G{uj,T) ^ g( 



T T* 



(25) 



with the Fermi liquid scale T* approximately given by 

T* ^TkKI ^ ( Ji - J2f , (26) 

where we introduced the renormalized anisotropy param- 
eter Kr as 



4(Ji~ J2) 
(Ji + J2)' 



(27) 



Throughout this paper we shall define T* as the energy 
scale at which Im t(w = T*) = 1.5 in the channel of larger 
coupling. Note that the parameter Kji can be consider- 
ably larger than the naive estimate, (Ji — J2)/{Ji + J2) 




K 



FIG. 5: Real part of the conductance computed from 
Eqs.(l22}, Eqs.(IT6j, and ([T7|l . The inset shows the universal 
collapse. 




(oX (~a/T^) 

FIG. 6: Imaginary part of the conductance from Eqs.i 
Eas. (|16|l . and p7p . The inset shows the universal scaling 
curve. 



due to the renormalization of the couplings Ji and J2 
in the high energy regime, D > oj > Tk- In the limit 
of T*, w < and T the conductance G{ui,T) be- 
comes a universal function of uj/T*, 



Gs 



-{uj,T^O)=GoF±iuj/T*) 



(28) 



The signs ± refer to the cases Ji > J2 and Ji < J2, 
respectively, and the scaling functions F±{y) have the 
properties 



a± + b± (^) 



2 

1/2 



l/2±c(^) , 



(29) 



In other words, for a; <C T* the conductance through the 
dot is Fermi liquid like, and Re G shows a ^ (uj/T*)'^ 
behavior, while for Tx ^ us ^ T* the real part of the 
conductance scales to its two-channel Kondo value with 
a small but increasing correction, ~ y/T* /uj. The latter 
behavior breaks down once the amplitude of the lead- 
ing irrelevant operator, ^ ^lu/Tk, reaches that of the 
anisotropy operator, ~ y^T* /uj, i.e. at frequencies in the 



yrV^, I.e. 
The constants a±, 6±, and c above 



range ui 

are numbers of order unity that depend somewhat on 
electron-hole symmetry breaking, but close to electron- 
hole symmetry a+ « 1, and a_ w 0. Note that the 
precise value of the constants b± and c depends also on 
the definition of the scale T*. 

The imaginary part of F±(y) has somewhat different 
properties and behaves as 




for w < T* 
for w > T* 



(30) 



In other words, the imaginary part of G must show a 
bump of size ~ Go at frequencies w ~ T*. These expec- 
tations shall indeed be met by our numerical results. 
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Similar to channel asymmetry, an external magnetic 
field also suppresses the non- Fermi liquid behavior, and 
introduces a new Fermi liquid scale, 



K 



(31) 



However, the magnetic field does not result in such a 
dramatic change in the conductance as the channel- 
symmetry breaking: while at w = the conductance 
exhibits a jump as a function of the channel- anisotropy, 
it changes continuously as a function of the magnetic field 
and shows only a cuspr^'^^ 

G{B)j,=j, « ^(l _ ;3 M \n{TK/\B\)) , (32) 

as it obviously follows from the singular behavior of the 
conduction electron phase shifts at the Fermi energy i^^i^^ 
As we shall see later, the AC conductance displays much 
more interesting features in a finite magnetic field. 



IV. NUMERICAL RESULTS 

In this section, we shall use the numerical renormaliza- 
tion group (NRG) method^^ to compute the eigenvalue 
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FIG. 7: (a) Imaginary part of the eigenvalue of the on-shell 
T-matrix, as function of lu/Tk, for several diflerent values of 
the anisotropy parameter, Kn = 4(Ji — J2)/(Ji -I- J2)^. In all 
cases Ji + J2 = 0.2. Curves with Ji > J2 or Ji < J2 scale 
to Im t{0) = 2 or Im t{0) — 0, respectively. The critical curve 
corresponding to Ji = J2 separates these two sets of curves, 
(b) Im t{Lj) for Ji = J2, as a function of yJuTfTK. The dashed 
line is a guide to the eye. (c) T* as the function of K\. 



of the T-matrix and from that the AC conductance. Al- 
though Eqs. (Uni) and (fTT]) hold at any temperature, finite 
temperature calculations are extremely delicate close to 
a non-Fermi liquid state. Therefore, we shall present nu- 
merical results only for T — Q temperature here. Nev- 
ertheless, according to the basic principles of scaling, a 
finite frequency uj plays a role rather similar to that of 
a finite temperature, and therefore the T = tempera- 
ture AC conductance, G{lo,T = 0), behaves rather sim- 
ilarly to the DC conductance at a finite temperature T, 
G(w = 0,T). 

To perform accurate calculations we assumed an 
electron-hole symmetrical conduction band and strongly 
exploited the symmetries of the Hamiltonian. The nu- 
merical results presented here have been obtained us- 
ing a new "flexible" NRG code, that handles symmetries 
dynamically 1^ In particular, in the absence of an external 
magnetic field, we used a symmetry S'C/ci(2) ® S'f7c2(2) ® 
SUs{2), with SUci{2) and SUc2i2) the charge SU{2) sym- 
metries in channels 1 and 2, respectively^^ and SUs{2) 
the spin SU (2) symmetry. The advantage of this sym- 
metry is that it is not violated even for Ji ^ J2, and it 
breaks down only to SUci{2) (g) SUc2{2) ® Usil) in the 
presence of a magnetic field. For the channel anisotropic 
cases we have retained a maximum of 750 multiplets dur- 
ing the NRG calculations, whereas 850 multiplets were 
kept in the presence of a magnetic field. All calculations 
were carried out with a discretization parameter A = 2. 
To compute the AC conductance, we have determined 
the composite Fermion's spectral function which, apart 
from an overall normalization factor, is equal to Im t(uj). 
This normalization factor can be easily fixed for Ji — Ji 
using the condition, Eq. (|20p . This procedure is much 
more accurate than estimating the normalization factor 
from the bare couplings, since the latter procedure suffers 
from the NRG discretization problem as well as from the 
loss of spectral weight at high energies, leading generally 
to a few percent error in the amplitude. 



A. Channel symmetry breaking 

First we investigated numerically how the non-Fermi 
liquid structure appears in the AC conductance through 
the double dot and how channel anisotropy destroys 
this non-Fermi liquid behavior. Some typical results are 
shown in Fig. [TJ for J\ ~ J2 we recover the two-channel 
Kondo result, Im t{uj ^ Q) = 1, and the deviation from 
the fixed point value scales as ~ -y/w / Tk , in agreement 
with Eq. 

For Ji 7^ J2 the new cross-over scale T* appears below 
which Im t(uj) crosses over from the two-channel Kondo 
value Im t(ti;) = 1, to lmt{u>) = 2 for Ji > J2 or to 
Im t(u;) = for Ji < J2 in the electron-hole symmetri- 
cal situation studied numerically. In the limit T* <^ Tk 
this cross-over is described by universal cross-over func- 
tions, similar to Eq. We determined these scaling 
functions numerically and displayed them in Fig. [5] (The 
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FIG. 8: Imaginary part of the on-shell T-matrix in the pres- 
ence of channel anisotropy as the function of oj/T* . The upper 
part corresponds to Ji > J2 while the lower part to Ji < J2. 
In both cases for T*,u! <^ Tk the curves follow the universal 
cross-over function, corresponding to a (aj/r*)^-like scaling 
at low frequencies and a 1 ± c {T* /uj)^^^ behavior at large 
frequencies. 



black curves were obtained by taking an extremely small 
value of T*, and chopping off the parts near lu ^ Tk-) 
The Fermi liquid scale T* extracted from t{uj) is shown 
in Fig. [71(c), and is in excellent agreement with the ana- 
lytical expression, Eq. (j^ . 

According to Eq. the real part of the conduc- 

tance can be computed from Im t{uj) through a simple 
integration. The resulting conductance curves are shown 
in Fig. [21 The behavior of Re G{uj) is strikingly simi- 
lar to that of Im t: it also exhibits a ^ ^/uJ singularity 
for Ji = J2 and crosses over from a value G = Go/2 to 
G = Go or to G = at the scale T* following the uni- 
versal cross-over functions, F±{lo /T*). We remark here 
that there seems to be no other reliable way than NRG 
to determine these universal cross-over functions, which 
connect two separate strong coupling fixed points, the 
non-Fermi liquid liquid fixed point and the Fermi liquid 
fixed point. These universal cross-over functions consti- 
tute some of the central results of this work. 

Performing a Hilbert transform, we also determined 
numerically the real part of the T-matrix, 'RBtiio), and 
from that the imaginary part of the conductance. These 
results are shown in Fig. [TOl It is quite remarkable that, 
although the scaling is not perfect because of the insuf- 
ficient accuracy of the Hilbert transform and the various 
integrations, clearly, the amplitude of the low temper- 
ature peak aX Lo ^ T* does not change as T* goes to 
0. (Note that T* varies over two orders of magnitudes.) 
This behavior is indeed expected based upon Eq. 
The numerical results confirm that for Ji > J2 and 
Ji < J2 the coefficients d± have different signs, d+ > 0, 
and d_ < 0, and that Im G(w) has a double peak struc- 
ture: it has one peak at w ~ Tk corresponding to the 
cross-over to the two-channel Kondo fixed point, and also 
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FIG. 9: (a) AC conductance as the function of uj/Tk- For 
Ji > J2 and Ji < J2 the curves scale Re G ^ Go and 
Re G — > 0, respectiv ely. In set: AC conductance for Ji = J2 
as the function of \JujITk- (b) AC conductance for positive 
(upper part) and negative (lower part) channel anisotropy pa- 
rameters as the function oiuj/T* . For u, T* <^ Tk, the curves 
follow the universal cross-over curves. 



another peak at lu T* related to the non-Fermi liquid 
Fermi liquid cross-over. 

It is interesting to observe from Figs. [51 [SI and [TUI that 
the range of two-channel Kondo scaling does not reach 
from Tk down to the cross-over scale T*, but rather it 
stops at a much higher energy scale, ~ ^/T*Tk, where 
corrections from the leading relevant operators start to 
dominate over the leading irrelevant operator of the two- 
channel Kondo fixed point. 



B. Effects of magnetic field 

We also performed calculations for Ji = J2 in the pres- 
ence of a local magnetic field, B. As mentioned earlier, a 
small local magnetic field destroys the non-Fermi liquid 
state and drives the system to a trivial, Fermi liquid fixed 
point below a scale Tb = B'^ /Tk < Tk^ 

Some typical results are shown in Fig. [Til At large 
magnetic fields, B > Tk, the Kondo resonance is clearly 
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FIG. 10: Imaginary part of the AC conductance as the func- 
tion of lu/Tk- Lower panel: Same as a function of lo/T* . 



split below the Zeeman field, and Re G{lo) exhibits a dip 
for |cij| < B. The width of this dip gradually decreases 
as one decreases the size of the field B, and its depth 
becomes smaller and smaller. However, it is not clear 
from the numerics if there is a critical field value. Be 
below which the dip actually disappears, as is the case, 
e.g. for the single-channel Kondo model. In fact, the nu- 
merical results seem to show just the opposite, i.e. that 
Re G{uj) remains a non-monotonous function in any fi- 
nite magnetic field, and only the height and width of the 
dip at w ~ Tb get smaller and smaller for smaller mag- 
netic fields while the dip itself is always present. This 
would indeed naturally follow from a simple scaling ar- 
gument: for B < Tk a magnetic energy scale is gen- 
erated, Tb = B^/Tk, and at this energy the real part 
of the conductance is expected to be Re G{lu « Tb) ~ 
Go [1/2 - a\B\/TK]. On the other hand, from Bethe 
Ansata^i we know the exact phase shifts, and from that it 
immediately follows that the DC conductance is given by 
Giij = 0) « Go [1/2 - C \B\/TK\og{TK/\B\)] at T = 
temperature, with G a constant of the order of unity3^ 
This observation suggests that in any finite magnetic field 
G(ijj) displays a dip, which has a width Auj ~ Tb, and 
height AG - \B\/TK\og{TK/\B\). Similar behavior is 
expected as a function of temperature, too. 

It is not clear either, if G{uj) becomes a universal func- 
tion of uj/Tb- In fact, it has been shown in a special, 
very anisotropic limit that no such universal function ex- 
ists for the non-linear DC conductance)^ We can argue 
that the same probably holds for the linear AC conduc- 
tance, although we do not have a rigorous proof. 

Unfortunately, from a numerical point of view the cal- 
culations in a magnetic field turned out to be extremely 
difficult: first of all, for the positive and negative fre- 
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FIG. 11: Top: Imaginary part of the on-shell T-matrix in 
the presence of a magnetic field and no channel asymmetry 
as the function ui/Tk- Lower curves correspond to higher 
magnetic fields. Bottom: AC conductance in the presence of 
a magnetic field and no channel asymmetry as the function 
of lo/Tk. Lower curves correspond to higher magnetic field 
values. 



quency parts of the spectral function one looses some- 
what different amounts of spectral weight. This effect 
turns out to be extremely large in the 2-channel Kondo 
case, probably as a consequence of the extreme sensitiv- 
ity of the non-Fermi liquid fixed point to the magnetic 
field. Therefore, for a given spin direction, one needs 
to match these positive and negative-frequency parts at 
the origin. Although this is a standard procedure fol- 
lowed by most groups, this leads to a large uncertainty 
in case of the 2-channel Kondo model. In fact, despite 
the extensive symmetries used, we were not able to ob- 
tain data of sufficient accuracy in the most interesting 
regime, lu <S^Tb = B^/Tk <C Tk, even using Hofstetter's 
density matrix NRG (DMNRG) method,^^ Therefore, we 
were not able to investigate the issue of universal cross- 
over functions for Ji = J2 and Tb = B'^/Tk < Tr-. We 
therefore consider these numerical results only as indica- 
tive but not decisive. 

We also need to recall the well-known fact that NRG 
produces an artificial broadening proportional to uj of 
the peaks occuring at finite frequencies. Thus, the cor- 
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rect shape of these spht peaks is presumably significantly 
sharper than that shown by the NRG results. 



V. CONCLUSIONS 

In this paper, we have studied the AC transport prop- 
erties of a double dot device realized recently by Potok 
et al. to reach the two-channel Kondo fixed point. First 
we derived an analytical expression for the linear conduc- 
tance in the regime where charge fluctuations are small 
and the system can be described by a Kondo Hamilto- 
nian. Our formula relates the AC conductance to the 
eigenvalue tiuo) of the dimensionless on-shell T-matrix, 
and is valid at any temperature and for any frequency. 
Our expression is the analogue of the formula obtained 
recently by Sindel et al. for the Anderson modcP'^ and it 
carries over to most Kondo- type Hamiltonians. 

The general properties of the two-channel Kondo fixed 
point, known from conformal field theory, allowed us to 
make many quantitative and qualitative predictions for 
the AC conductance, G{uj): for equal couplings to the 
two channels, G{uj) shows an ^Juj/Tk singularity at the 
two-channel Kondo fixed point. Using the results of con- 
formal field theory,'^*' we were able to compute the real 
and imaginary parts of the function G{lu,T) and deter- 
mine the corresponding scaling functions for both the real 
and the imaginary parts of the conductance through the 
dot in the universal regime, oj, T <C Tk and Ji = J2. The 
generic properties of the AC conductance in this regime 
are summarized in Fig. 1121 

Conformal field theory also gave us a way to predict the 
basic properties of Re G{lu) and Im G{uj) at T = tem- 
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FIG. 12: Sketch of the real and imaginary parts of the AC 
conductance for Ji = J 2 and u>,T <^ Tk- 



FIG. 13: Sketch of the real and imaginary parts of the T = 
temperature AC conductance for Ji > J2. The various powers 
shown in the picture follow from conformal field theory. The 
high frequency behavior is a result of perturbation theory. We 
assumed electron-hole symmetry. 



perature, in the presence of channel symmetry-breaking 
(see Fig. [131). For Ji ^ J2 Re G{uj) crosses over to a much 
smaller or a much larger value (depending on the sign of 
asymmetry) at the Fermi-liquid scale T*, below which it 
becomes an analytical function of uj. This cross-over at 
Lo T* is described by universal cross-over functions that 
we have determined numerically. The asymptotic prop- 
erties of the real and imaginary parts of the conductance 
are dictated by conformal field theory (see Eqs. (|29p and 
([50]) ). It is quite remarkable that Im G{ijj) shows a dou- 
ble peak structure at frequencies lo ^ T* and uj ^ Tk- 
Both peaks are of amplitude ~ Go but the sign of the 
peak at a; '--^ r* changes with the sign of Ji — J2- 

One of the important conclusions that one could draw 
from the analysis of G{uj) was, that the two-channel 
Kondo regime is, in a sense, restricted to the regime, 
y/TxT* <T,uj <Tk-- Although it is true that the entire 
regime, T* < T,uj < Tk is governed by the 2-channel 
Kondo fixed point, for T,ui < \/TkT* the leading rele- 
vant operator is more important than the leading irrele- 
vant operator, and therefore, the scaling curves charac- 
teristic to the two-channel Kondo fixed point itself cannot 
be seen in this regime. This refines somewhat the phase 
diagram of the two-channel Kondo model, as already in- 
dicated in Fig. [21 The two-channel Kondo scaling regime 
is thus limited by a boundary ~ | Ji — J2 1 ■ 

We have also investigated the effects of a small Zee- 
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man field on the AC conductance. For B > Tk the AC 
conductance exhibits a dip whose width is just B. Nu- 
mericahy we find that, apparently, this dip survives for 
any small magnetic field, B < Tk- This would indeed 
be in agreement with a simple scaling argument we pre- 
sented, that also predicts a similar behavior as a func- 
tion of temperature. In other words, at the two-channel 
Kondo fixed point the Kondo resonance appears to be 
split at any magnetic field. Unfortunately, while our nu- 
merics seems to support this picture, it is not accurate 
enough in the regime, B <C Tr- to give a decisive answer. 
We remark that the logarithmic magnetic field depen- 
dence of the phase shift would also probably imply that 



universal scaling (i.e. T/Tb scaling and the disappear- 
ance of the scale Tk for T, Tb ^ Tk) should be destroyed 
by logarithmic corrections in the presence of magnetic 
field. 
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